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Production of J/ip from nucleus- nucleus reactions de- 
pends sensitively on the dissociation cross section with light 
hadrons. Effective lagrangian methods are used to describe 
the hadronic degrees of freedom, including strangeness and 
charm. Cross sections with pions, rho mesons, kaons and nu- 
cleons having magnitudes 4-8 mb are found, and with steep 
thresholds. This, folded with thermal momentum distribu- 
tions for the scattering partners, suggests a mean dissocia- 
tion lifetime ~ 20 fm/c. Therefore, the "abnormal" J /ip sup- 
pression seen in recent Pb+Pb experiments seems to owe to 
expected hadron kinetics. 



Response of nuclear matter to high energy densities 
affords the possibility of creating in the laboratory a sys- 
tem of quantum chromodynamic (QCD) matter, the so- 
called quark gluon plasma. Among the signatures for 
creation of the mesoscopic colored volume is an idea put 
forward in 1986 by Matsui and Satz M to look at elec- 
tromagnetic spectra for evidence of charmonium, bound 
states of charm-anticharm. They are very tightly bound 
and consequently relatively small hadrons which ought 
to effectively probe superdense hadronic matter. Their 
utility in this context comes from the fact that Debye 
[color] screening in the plasma would so strongly sup- 
press cc binding in favor of charm propagating decou- 
pled from anticharm to later join with more abundant 
light anti-quark species forming D mesons, that suppres- 
sion of J /ip would indicate plasma formation. Above 
some critical temperature, the screening radius becomes 
smaller than the binding radius offering this possibility. 
J /ip suppression has since been regarded as a promising 
signature of quark gluon plasma formation. 

Charmonium production cross sections from proton- 
induced reactions that are directly proportional to the 
target mass number A, are said to be normal. And 
yet, proton- nucleus and nucleus-nucleus experiments per- 
formed over the past several years have exhibited a com- 
mon A a dependence, with a « 0.92 |-|]. This suppres- 
sion is now understood as being due to absorption of the 
precursor state to J /ip on nucleons and so is in some 
sense of normal hadronic consequences. Additional sup- 
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pression of ip' production revealed in nucleus-nucleus ex- 
periments has been attributed to absorption on produced 
comoving hadrons |^|. However, recent experimental 
studies at CERN have investigated J/ip production in 
Pb+Pb reactions at 158 GeV/nucleon and have reported 
a dramatic "abnormal" suppression compared to lighter 
systems' (AB) a behavior 0. It is reasonable to begin 
discussing the possibility of deconfmement being respon- 
sible. But first, all conventional suppression mechanisms 
must be under control. 

Absorption of J /if) or ip' on comovers has been pro- 
posed as the dominant dissociating effect. One mecha- 
nism owes to deconfmement; if the plasma were present, 
charmonium states would not form at all due to screen- 
ing. This is simply a restatement of the signature idea of 
Matsui and Satz. Another mechanism discussed in the 
literature is pre-thermal dissociation, where the charmo- 
nium is excited above the (DD) threshold either by a 
partonic medium which has not had sufficient time to 
equilibrate |8|, or by color flux tubes in the infant QCD 
medium [p|flC|] . The effectiveness of pre-thermal suppres- 
sion mechanisms rests on separation of time scales: the 
first mechanism requires that dissociation proceeds faster 
than thermalization and the second requires that color 
flux tubes decay more slowly than dissociation. Finally, 
we come to thermal dissociation where the (cc) is ab- 
sorbed through one of the myriad processes involving a 
light hadron. 

The crucial question is the magnitudes of cross sec- 
tions for J /ip + h, where h is one of the set (it, p, K, 
N . . .). Estimates up to now have included 1) a pertur- 
bative approach at the quark level where light hadrons 
interact with J /ip solely through gluonic content effects 
pUP^I , 2) a nonperturbative study with quark exchange 
including a confining potential |l3f| , and finally 3) an ef- 
fective lagrangian approach Jli} |. When the perturbative 
approach of Peskin and Bhanot was applied by Kharzeev 
and Satz to pions interacting with J /ip, a cross section 

~ 0.1 mb was found for y/s roughly 1 GeV above thresh- 
old p5|| . Results for the nonperturbative calculation of 
Martins, Blaschke and Quack applied to pions revealed 
cross sections peaking at several mb also about 1 GeV 
above threshold. The hadronic field theory calculation 
of Matinyan and Miiller resulted in cross sections with 
pions and rho mesons of order 1 mb for similar energies. 
This situation is unsettling as there is roughly two or- 
ders of magnitude discrepancy in these estimates. The 
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aim of this letter is to report on further investigation of 
the issue using effective lagrangian methods quantifying 
cross sections within a consistent gauge invariant treat- 
ment (including contact terms which are missing in pre- 
vious analyses) and including a more complete hadronic 
medium as input to kinetic theory for a baseline estimate 
of the dissociation rate in hadronic matter. 

Hot hadronic matter is populated most abundantly 
by 7r, K, and p If circumstances of an apprecia- 

ble baryon chemical potential arise, the nucleons become 
important as well. Each species can induce charmonium 
dissociation with respective final states governed by con- 
servation laws. A consistent treatment of light mesons, 
heavy mesons and possibly baryons is needed to work to- 
ward a reasonably reliable description for hadronic mat- 
ter. We begin with an SU(4) symmetry so as to include 
charm, and we introduce pseudoscalar (P = if a X a ) and 
vector ( — v£\ a ) meson matrices, where ip a and v% 
are pseudoscalar and vector multiplets and the As are 
SU(A) generators [0. The symmetry is severely broken 
due to the large charm quark mass, so we use physi- 
cal mass eigenstates and matrices, and incorporate con- 
straints where possible. 

The free meson lagrangian reads 

£ = Tr(<9' i P t «9 A( P) 



+mass terms. 



(1) 



Interactions are generated by replacing the space- 
time derivative with a gauge covariant one d^f — > 
DJ=dJ+ [A^f], where A„ = -ig/2V„. The vector 
mesons are recognizable as playing roles of gauge bosons. 
As is usual in effective field theory strategies, to keep 
gauge consistency, we must collect terms of order g 2 . 
They are (using = P and = V) 

= ijTrfPFa/ - d^PV^P) 

+ -g 2 Tr(PV^ V^P-PV^PV^) 

+ igTY{d» V V [V^ V V ] + [W, V v ]d^V v ) 
+ g 2 Tr(V' t V \V V„]). (2) 

The order g terms deliver three-point vertices and the or- 
der g 2 terms are responsible for so-called contact terms, 
or four-point couplings which are necessary for a gauge 
invariant theory. The interactions of interest for this 
study are the following (using ifj as shorthand for J/ip) 

i^DD* = \g-*DD* (DTiD^d^i-d^DnDliti 
-H.c), 

% — 

£ p d*d* = -^g P D>D*P? (D^nidf.Dl) - (d^D* u )nDl 



-(d u D;) Ti D* u + (dvDDnD*") , 

Udd = ig^DDV (Dd^D - (d M D)D) , 

-(d v D;)D* v + {d v Dl)D*»), (3) 

and 

^4>ttDD* — —g^pDD g-KDD'^ D* ^TiDm, 
^t/jpDD = -g P DD g-4,DDlp^ Pn,iDnD, 

i-^pD'D* = -g^D-D* g P D*D*if> tl (2p flt iD* I/ TiD* v 



(4) 



The K-D s -D* interaction has the same structure as the 
ir-D-D*, complete with contact terms. 

In this approach there are several coupling constants 
in Eqs. (0) and (|]). Methods for constraining them will 
be discussed below. The 0{g 2 ) terms carry one power 
of coupling constant for each three-point vertices from 
which the contact term collapses. As a specific example, 
the direct, exchange, and contact graphs for reaction J/ip 
+ it — > D* + D are shown in Fig. [t] (a), (b) and (c), 
respectively. Contribution to the amplitude from each 
graph is proportional to g^ DD ,g^ DD . 
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FIG. 1. Feynman graphs for the process J /t[> + n — > D + 



D* 



Appealing to measured hadronic decays to constrain 
the coupling constants gains very little in this case as the 
best measurement of the D*^ width has an upper limit 
of 131 keV jL8|. Model calculations based on relativistic 
potential approaches suggest T D ,± =46 keV This 
corresponds to a value for g^DD* of 8.8, which will be 
used here. In the absence of empirical constraints, all 
other coupling constants are obtained from vector meson 
dominance or flavor symmetry arguments. Vector domi- 
nance gives Q 



g P DD = g P d*d* 



5.6, 



9i/>DD — 9ipD*D* 



7.7. 



(5) 
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For the kaons, we use SU (4) symmetry to relate gx d s d* 
— V^gpTTTT = 8.5, and the pirn coupling is adjusted to 
give a p width of 151 MeV. 

The following absorption reactions are considered 
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d + d* s , d + d; 



(6) 
(7) 

(8) 

Just as in Fig. [l], for a reaction of type (^), each of the 
reactions listed in (||)-(||) has a direct, exchange and a 
"seagull" graph contributing to the amplitude. Full ex- 
pressions for the amplitudes and other details will be pub- 
lished elsewhere EjJ. 

Results for cross sections as functions of y/s are pre- 
sented in Fig. ||. A startling feature in the endother- 
mic reactions is the sharp rise just above threshold. The 
cross sections reach maximum strength a few hundred 
MeV above their respective thresholds. We note the sig- 
nificant difference between these results and the ir- and 
p-induced reactions from Ref. |l4[ , which are due to the 
inclusion of contact terms and interference effects. Fur- 
ther, we note the relatively large kaon cross section. 
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FIG. 2. Dissociation cross sections for tt, p, and kaon. Nu- 
cleon-induced dissociation is also included. See text for de- 
tails. 



Next we turn to the baryons and begin by writing an 
effective lagrangian describing the interaction with pseu- 
doscalars. The full SU (4) expression is quite lengthy , 
so we concentrate on the particular term allowing analy- 
sis of the reaction J/tp + N^A c + D. The lagrangian 
bearing pseudoscalar coupling is p2j 



The coupling constant is constrained by QCD sum rules 
which predicts a value gDNA c = 6.7 ±2.1 There are 
two Feynman graphs, a direct and exchange contribution, 
which when added together once again produces an am- 
plitude respecting gauge invariance. The cross section is 
presented also in Fig. |[ It to rises sharply at threshold 
to ps 7 mb and falls with rising energy. 

After the early stages of a high-energy heavy-ion reac- 
tion, the produced hot and dense system rapidly expands 
and cools leaving a kinetically thermal and probably even 
chemically equilibrated hadronic fireball. The fireball 
quickly cools, most likely falls out of chemical equilib- 
rium |2jj2^], until it eventually freezes out. Therefore, 
charmonium kinetics over a range of temperatures are 
needed for establishing predictions for J/ip production. 
Relativistic kinetic theory allows for simple estimates of 
rates. Nonequilibrium effects are of course important, 
and will be discussed elsewhere 12(1. In general, the in- 
variant four rate for a + b scattering assuming Boltz- 
mann distributions is 



dR 



ill 
(2n) 



dz\(z 2 T 2 ,m 2 a ,m 2 b )daK 1 (z), (10) 



where g is an overall degeneracy factor, z = y/s/T, z m i n 
= max(m a +mb,M) / T, with M being the sum of all fi- 
nal state particles' masses, and K\ is a modified Bessel 
function. The number of times particle a scatters with a 
particle of type b per unit time is then 



Rate = 



dR 
dN~ b : 



where dNf, is the number density of b particles, 
dN b = -^-Tm 2 K 2 (m b /T). 



(11) 



(12) 



Again, K2 is a modified Bessel function. 

The number of J /ip dissociations per unit time induced 
by each light hadronic species are shown separately in 
Fig. |^ as well as a combined sum. J ftp dissociation rate 
in resonance matter is w 0.03 (fm/c) -1 . If we look to- 
ward temperatures approaching 200 MeV, the rate ap- 
proaches 0.1 (fm/c) -1 . It was previously thought that 
thermal hadronic dissociation rates were so small as to 
be insignificant on the time scale of the fireball created in 
heavy ion reactions. Fig. ^ indicates that reactions with 
light hadrons will indeed be an important hindrance for 
charmonium production. 



Hdna c = g DN A c (N l5 A c D + A cl5 ND) . 



(9) 
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FIG. 3. The thermal rates for J /ip dissociation induced 
by 7r (short-dashed curve), p mesons, (dotted curve), kaons 
(dot-dashed curve), and nucleons (long-dashed curve), and 
the sum (solid curve). 



There are of course error bars to discuss associated 
with these calculations. First, the coupling constants are 
uncertain to a level of ten of percent or so. Form factors, 
which would tend to reduce the cross sections, have not 
been included. On the other hand, a long list of reactions 
involving other mesons could be imagined, and some of 
which might be important. For instance, axial vector 
charm mesons Di(2420) |26| have relatively large widths 
and could therefore play an important role in dissociating 
J through such reactions as 



7T + J/V>-> D 1 +D\ D 1 + D* 



(13) 



Another candidate likely affecting J/ ip dissociation in the 
medium is Xc- 

Finally, the rate estimates made here would be in- 
creased dramatically if phase space were overpopulated, 
interpretable as finite chemical potentials. As the fireball 
expands isentropically, pion and kaon chemical potentials 
of order 50-100 MeV develop in model calculations po] . 
This would easily buy back a significant factor in the 
rates. All these effects are currently under investigation 
and will be reported upon separately p3]. 
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